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Abstract. Quantum Chromodynamics (QCD) is generally assumed to be the fundamental 
theory underlying nuclear physics. In recent years there is progress towards investigating the 
nucleon structure from first principles of QCD. Although this structure is best revealed in Deep 
Inelastic Scattering, a consistent analysis has to be performed in a fully non-perturbative scheme. 
The only known method for this purpose are lattice simulations. We first sketch the ideas of 
Monte Carlo simulations in lattice gauge theory. Then we comment in particular on the issues 
of chiral symmetry and operator mixing. Finally we present our results for the Bjorken variable 
of a single quark, and for the second Nachtmann moment of the nucleon structure functions. 



1. Introduction 

For 37 years we have a strong candidate for a fundamental theory underlying nuclear physics. 
There is a well-founded paradigm that the whole of nuclear physics can in principle be explained 
by QCD. Nevertheless this theory has hardly been represented at this workshop, because most 
talks referred to complex phenomena, which can in practice still not be computed based on 
QCD. Effective approaches — like the shell-model — are popular for this purpose. They are 
often successful for specific questions, but they all have their limitations. 

In this contribution we do address QCD, and we report in particular on progress in the 
understanding of the nucleon structure function from first principles. This is obviously a non- 
perturbative task, and the only method to tackle it systematically are lattice simulations. 

The formulation of QCD looks amazingly simple: it deals with quarks with a local SU (3) 
symmetry ("colour symmetry"), which gives rise to strong interaction mediated by 8 gauge 
bosons ("gluons"); one for each generator of the gauge group SU(3). Since the latter is non- 
Abelian, the gluons also interact among each other, which makes the dynamics of the system 
complicated. In the low energy regime, which dominates our daily experience, the gauge coupling 
is strong, and the interactions are extremely involved. This can be seen for instance from the 



fact that the fundamental masses of the lightest (valence) quark flavours — which we assume 
to be generated by the Higgs mechanism — only contribute about 2 % to the nucleon mass; all 
the rest is contained in some dense "foam" of gluons and sea-quarks. 

So far perturbation theory could be carried out for some quantities up to af , with a s = gl/Air, 
g s being to strong gauge coupling. This is applicable to a number of high energy processes, 
dominated by the exchange of few gluons. A suitable method for that purpose is dimensional 
regularisation, which is, however, strictly limited to perturbation theory. 

On the other hand, the lattice regularisation enables a fully non-perturbative treatment, i.e. 
we can really address finite values of a s . Here we refer to the functional integral formulation 
of quantum field theory, which is (in principle) based on an integration over all possible 
configurations of the fields involved. The partition function, and all n-point functions that can 
be derived from it, are then given by an integration over all these configurations with a phase 
factor exp(iS'), where S is the action of the configuration. The convergence of such integrals 
can be improved drastically by moving to Euclidean space-time, i.e. by invoking an "imaginary 
time" (Wick rotation). This changes the factor mentioned above to exp(— S'e), where S'e is 
the Euclidean action. This transition is justified in view of the expectation values of n-point 
functions if four conditions — known as the "Osterwalder-Schrader axioms" — hold. For many 
theories, including QCD at zero vacuum angle 9 and zero baryon density, S'e is real positive for 
all configurations, hence the factor exp(— Se) can be interpreted as a probability distribution 
for the possible field configurations. 

Lattice studies consider discrete points in a finite volume of Euclidean space-time. The 
matter fields — in this case the quark fields \& — are located on these lattice sites, while 
the gauge field variables E SU(3) live on the links connecting nearest neighbour sites. Thus 
the functional integral is completely well-defined. Most theories of interest are bilinear in the 
fermion fields. Then the action takes the structure Se^, U] = ^M[U]^+S ga _ uge [U], where the 
indices of ^, ^ run over all lattice sites, and on each site over all internal degrees of freedom. 
The fermion matrix M contains the gauge couplings between neighbouring sites (due to the 
discrete covariant derivative in the Dirac operator). S ga _ uge is the pure gauge part of the action. 
In accordance with the Spin-Statistics Theorem, the components of the spinor fields Vl/ are 
given by Grassmann variables, i.e. they anti-commute, and the corresponding integration yields 

Z = J V^V^VU exp(-5 E [*,*,tf]) = J VUdetM[U]e- s ^ usc[u] . (1) 

The functional measure T>U now represents the integral over the (compact, gauge invariant) 
Haar measure for each link variable. Thus we preserve exact gauge invariance on the regularised 
level. For system sizes of interest, it is still far beyond our ability to compute this integral 
explicitly. Hence the method of choice is to generate a large set of gauge configurations with 
the given probability distribution; moreover the configurations should be well de-correlated, i.e. 
statistically independent from each other. There are powerful algorithms for this purpose, see 
Refs. pQ for reviews. Once such a set of configurations is generated, it can be used to measure 
observables of interest. This approximate functional integral yields results for expectation values 
with statistical errors (due to the finite sample of configurations) and systematic errors (finite 
lattice spacing, finite volume etc.). These errors can be estimated and reduced with additional 
computational effort. The important point, however, is that the result is fully non-perturbative. 
This allows us nowadays for instance to calculate the light hadron spectrum from first principles 
of QCD up to a precision of a few percent. The results agree with phenomenology within these 
errors. 

Our daily experience is dominated by nuclear physics at low energy. Regarding QCD this 
implies a dominance by the light quark flavours, in particular the flavours u and d. Their masses 



are very small compared to Aqcd — the intrinsic scale of of QCD, which provides a scale for the 
hadron masses. Hence these quarks have an approximate chiral symmetry, which is of prominent 
importance in low energy QCD. To briefly review this property, we first note that fermion fields 
can be decomposed into "left- and right-handed" components by means of the chiral projectors, 

* LjR = i(l±75)* , *L,R = t^(lT75) • (2) 

In the chiral limit (zero quark masses) these two components decouple. It is allowed to add 
quark masses to the Lagrangian, since QCD is a vector theory — the left- and right-handed 
components couple in the same way to the gauge field — (unlike the electroweak sector of 
the Standard Model). From a perturbative perspective, the approximate chiral symmetry still 
protects the masses of light quarks from large renormalisation. Based on this property, it is 
occasionally claimed that there is no hierarchy problem for fermion masses. However, it is the 
non-perturbative level that ultimately matters for physics, and there it is highly non-trivial to 
implement (approximate) chiral symmetry in a regularised system. A quite obvious ansatz is 
the Wilson lattice fermion, which starts from the naive discretisation of the derivative in the 
Dirac operator, and subtracts a discrete Laplacian to send the fermion doublers to the cutoff 
scale. However, this extra term breaks chiral symmetry explicitly. As a consequence, under 
gauge interaction the fermion masses are renormalised such that they naturally end up at the 
cutoff scale as well. The (observed) property m Uj d <C Aqcd can then only be attained by a 
tedious fine-tuning of negative bare quark masses, i.e. the hierarchy problem is back. 

Only at the very end of the last century was there a breakthrough in the formulation of chiral 
lattice fermions, which overcomes this problem, at least in vector theories. A key observation 
was that a lattice modified version of chiral symmetry of the form [2] 
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is sufficient to rule out additive mass renormalisation. Here a is the lattice spacing, and in the 
continuum limit this turns into the standard chiral transformation (infinitesimal form). Even 
on the lattice the fermionic Lagrangian ^fD^f is invariant, if the lattice Dirac operator D fulfils 
the Ginsparg- Wilson relation (GWR) 

{D, 75} = D l5 + 75 L> = aD l5 D . (4) 

On the other hand, the functional measure T>^T>^f is not invariant, and this is just what it takes 
to reproduce the axial anomaly correctly [2]. 

Solutions to the GWR can be constructed based on block spinor renormalisation group 
transformations ("perfect lattice fermion") OS], or by pulling apart the chiral modes in an 
extra "dimension" ("Domain Wall Fermion") [5j. By integrating out that "dimension", one 
obtains — for 75-Hermitian operators, = 75.D75 — the "overlap Dirac operator" D ov [6], 



aD ov = 1 + A ov , A ov = A/VaTA , A = aD - 1 . (5) 

It is easy to verify that it satisfies the GWR, which can now be written as D ov + Dl v = aDl v D ov . 
For the kernel D we can insert for instance the Wilson operator [6], but properties like scaling, 
locality and approximate rotation symmetry can be improved by using a better kernel [7]. 

In addition to the exact, lattice modified chiral symmetry, the GWR also prevents uncon- 
trolled operator mixing on the regularised level This powerful property is very helpful, for 
instance for the numerical study of a fully non-perturbative Operator Product Expansion (OPE). 



We now summarise our lattice OPE project [HI [10] as an example for the application of 
chiral quarks in lattice QCD. This takes us beyond the hadron spectrum; we investigate the 
internal structure of the nucleon (see Ref. [11] for a recent review). Phenomenologically this 
structure is revealed in Deep Inelastic Scattering of a lepton on a nucleon target. This scattering 
is dominated by a one-photon exchange involving a single quark, as sketched below. 
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Although this is a high energy process, there are problems with a perturbative treatment 
due to ambiguous UV and IR divergences. A fully satisfactory study based on QCD has to 
be consistently non-perturbative, i.e. it is a challenge for the lattice. The involvement of a 
single quark implies that chirality is conserved, hence the use of Ginsparg- Wilson quarks is 
appropriate. In our simulations we applied 2 flavours of degenerate overlap valence quarks with 
masses of 29 MeV and 73 MeV, which correspond to pion masses of 280 MeV and 440 MeV, 
respectively. The goal is the evaluation of the product of electromagnetic currents J M — of the 
photon exchange, with large transfer momentum q — between quark states of low momentum 
p. This product is decomposed by the OPE as follows, 

W, v ~ (9(p)\J^q)M-q)\9(p)) = £ ^ i ... Mn (g)(*(p)|og... M J^(p)) . (6) 

m,i,n 

are Wilson coefficients, which are independent of the target and therefore of p, and C>( m ) 
are local operators, which are relevant to describe the nucleon structure. They involve the quark 
momentum components p^.. We further use the Clifford index i = 1 ... 16, and the index m to 
distinguish operators with the same symmetries. 

A truncation of the OPE on the right-hand-side of eq. ([6]) requires the scale separation 

P 2 « q 2 • (7) 

We probed three photon momenta, \q^\ = 2.2 GeV, \q^\ = 3.3 GeV, \q^ \ = 4.4 GeV. Based 
on the above scale separation, we truncated the quark bilinears at 0(|p| 3 ). Thus we still include 
1360 operators, but for photon momenta q oc (1, 1, 1, 1) they occur in only 67 equivalence classes. 
The strategy is now to measure the elements of and (0( m ' ) ) for a multitude of low quark 
momenta p, and solve eq. ([6]) for the Wilson coefficients. They are over-determined in this way, 
but by means of a Singular Value Decomposition we arrived at reliable values [8]. In Fig. [I] we 
compare the results obtained numerically at photon momentum q^ to the Wilson coefficients at 
tree level in the chiral limit. We see that they follow the same pattern, though gauge interaction 
reduces the absolute values significantly. In particular the coefficients C\ and C-j . . . C\% vanish 
in the chiral limit, because they are attached to operators with an even power of 7- matrices. We 
obtain small absolute values for them also in the interacting case at finite quark mass, which 
confirms that chiral symmetry on the lattice works, and it suppresses operator mixing. 

As a first application, we evaluate the Bjorken variable x for one unpolarised quark. The 
expectation value (x) measures the nucleon momentum fraction carried by the quark involved 
in the Deep Inelastic Scattering. To be explicit, we consider the tensor [9] 

T„ v = -TrlS^^W^ip, q)] = S^Wx + p^p v W 2 + {p^q v + p v q^)W 4 + q^q v W 5 , (8) 
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Figure 1. The Wilson coefficients for quark masses m q = 0.011 and m q = 0.028 in lattice units, 
corresponding to 29 MeV and 73 MeV. In both cases the photon momentum amounts to 3.3 GeV. The 
results for the 67 independent coefficients are compared to their values on tree level. We see a similar 
pattern, with reduced absolute values under interaction. The coefficients which vanish in the chiral limit, 
C\ and CV . . . Cie, are still small at finite quark mass, which confirms the validity of our method. 

where S is the quark propagator. The Dirac-index structure is not visible here due to the trace 
(terms with W3 and We would occur for polarised targets). 

In particular T w is directly related to the Bjorken variable in the high energy Bjorken limit. 
In this limit, the Minkowski space momenta obey (p • q) 2 3> p 2 q 2 . Among the 22 operators in 
the expansion of T„„ we focus on those with spin 2, where the Bjorken limit leads to [9] 
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Fig. [2] (on the left) shows the numerator against the denominator of the right-hand-side in eq. 
([9]). A broad range of points measured at m q = 73 MeV is consistent with a linear behaviour, 
which allows us to read off (x) from the slope (the errors are statistical only), 
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Finally we proceed to the nucleon structure function |10| . On a lattice of spacing a, a generic 
moment can be expanded as 



M(q z ) = c 2 (q 2 a 2 )A 2 (a 2 ) + 



C4,{q 2 a 2y 



-A 4 (a 2 ) + ... 



(11) 



where A n are reduced nucleon matrix elements (the Lorentz structure is factored out), c n are 
the corresponding Wilson coefficients, and n is the twist. Both factors are bare quantities. If one 
tries to replace a by some scale \i of a soft renormalisation scheme, one runs into ambiguities, 
which do not occur on the lattice with its hard momentum cutoff ir/a. 

In our study, we evaluated the (spin averaged) nucleon matrix elements for all 67 types of 
operators, and extracted moments of the structure function by Nachtmann integration of the 
hadronic tensor [12]. This allows us in particular to project out operators of spin 2. For instance, 



the traceless form of W^ v leads to 
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Here n is a constant 4- vector with n 2 = 1, and the last expression in eq. (|12p refers again to the 
Bjorken limit. By means of a further projection [10] we can separate the structure functions F 2 




Figure 2. On the left: The term q 2 T^^ for spin 2, versus 2^ M>J ,p AI p„ (numerator and denominator 
in eq. ©). The squares, circles and triangles refer to photon momentum gW, q^ 2 \ q^\ respectively. 
In all three cases the data points are compatible with a linear behaviour. Thus the slope allows for the 
determination of the Bjorken variables (x) given eq. (llOp . On the right: The nucleon moment M.2 against 
the photon momentum squared. We distinguish the contributions of twist 2 and twist 4. In all cases 
the result for our two quark masses agree within statistical errors. For large q 2 the twist 4 contribution 
seems to vanish, whereas the twist 2 contribution stabilises at a finite value in the range w 0.5 . . . 0.8. 



and Fl- Based on the numerical measurement of the matrix elements, we can carry out the first 
integral in eq. (|12|) . Thus we arrive at explicit and non-perturbative results for A4. 

Our results for the second moment M2 are shown in Fig. [2] (on the right). They agree within 
the errors for both quark masses, at all three photon momenta, which suggests that they may 
be valid up to the physical pion mass. As we increase the photon momentum (i.e. the transfer 
momentum in the hard scattering), the twist 4 contribution tends to vanish, while the twist 2 
term stabilises at a moderate value. As usual the numerical errors increase for lighter quarks. 

These results demonstrate that the explicit determination of the nucleon structure functions 
- as observed in Deep Inelastic Scattering — directly from QCD is feasible, in the framework 
of a non-perturbative regularisation free of ambiguities. Moreover we obtained explicit results 
for the Bjorken variable of the quark involved in this scattering. Our results provide insight into 
the nucleon structure. For a detailed discussion we refer to Ref. [10] . 
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